COARSE-GRAINING SCHEMES FOR STOCHASTIC LATTICE 
SYSTEMS WITH SHORT AND LONG-RANGE INTERACTIONS 



MARKOS A. KATSOULAKIS*, PETR PLECHAC 1 ", LUC REY-BELLET*, AND DIMITRIOS 

K. TSAGKAROGIANNIS 5 

Abstract. We develop coarse-graining schemes for stochastic many-particle microscopic models 
with competing short- and long-range interactions on a ci-dimensional lattice. We focus on the coarse- 
graining of equilibrium Gibbs states and using cluster expansions we analyze the corresponding 
renormalization group map. We quantify the approximation properties of the coarse-grained terms 
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semi-analytical numerical schemes that are accompanied with a posteriori error estimates for coarse- 
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1. Introduction. Many- particle microscopic systems with combined short and 
long-range interactions are ubiquitous in a variety of physical and biochemical systems, 
[35] . They exhibit rich mesoscopic and macroscopic morphologies due to competition 
of attractive and repulsive interaction potentials. For example, mesoscale pattern 
formation via self-assembly arises in heteroepitaxy, [33], other notable examples in- 
clude polymeric systems, [TJ], and micromagnetic materials, [TB]. Simulations of such 
systems rely on molecular methods such as kinetic Monte Carlo (kMC) or Molecular 
Dynamics (MD). However, the presence of long-range interactions severely limits the 
spatio-temporal scales that can be simulated by such direct computational methods. 

On the other hand, an important class of computational tools used for accelerat- 
ing microscopic molecular simulations is the method of coarse-graining. By lumping 
together degrees of freedom into coarse-grained variables interacting with new, ef- 
fective potentials the complexity of the molecular system is reduced, thus yielding 
accelerated simulation methods capable of reaching mesoscopic length scales. Such 
methods have been developed for the study and simulation of crystal growth, surface 
processes and polymers, e.g., [19j [TT1 [25] Q] HI HI] , while there is an extensive litera- 
ture in soft matter and complex fluids, e.g., (39] [28l [HI H2] . Existing approaches can 
give unprecedented speed-up to molecular simulations and can work well in certain 
parameter regimes, for instance, at high temperatures or low densities of the systems. 
On the other hand important macroscopic properties may not be captured properly 
in many parameter regimes, e.g., the melt structures of polymers, [25] ; or the crystal- 
lization of complex fluids, [32] ■ Motivated in part by such observations we formulated 
and analyzed, from a numerical analysis and statistical mechanics perspective, coarse- 
grained variable selection and error quantification of coarse-grained approximations 
focusing on stochastic lattice systems with long-range interactions, [53] [5] [2T] . We 
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have shown that the ensuing schemes, known as coarse-grained Monte Carlo (CGMC) 
methods, perform remarkably well even though traditional Monte Carlo methods ex- 
perience a serious slow-down. In this paper we focus on lattice systems with both 
short and long-range interactions. Short-range interactions introduce strong corre- 
lations between coarse-grained variables and a radically different approach needs to 
be employed in order to carry out a systematic and accurate coarse-graining of such 
systems. 

The coarse-graining of microscopic systems is essentially a problem in approxima- 
tion theory and numerical analysis. However, the presence of stochastic fluctuations 
on one hand, and the extensive nature of the models (the presence of extensive quan- 
tities that scale as O(N) with the size of system N) on the other create a new set of 
challenges. Before we proceed with the main results of this paper we discuss all these 
issues in a general setting that applies to both on-lattice and off-lattice systems and 
present the mathematical and numerical framework of coarse-graining for equilibrium 
many-body systems. 

We denote by a microscopic states of a many-particle system and by Sn the set of 
all microscopic states (i.e., the configuration space). The energy of a configuration is 
given by the Hamiltonian Hn(<?) where N denotes the size of the microscopic system. 
An example studied in this paper is the <i-dimcnsional Ising-type model defined on a 
lattice with N = n d lattice points, and suitable boundary conditions, e.g., periodic. 
For both on-lattice or off-lattice particle systems the finite- volume equilibrium states 
of the system are given by the canonical Gibbs measure at the inverse temperature 
ft, describing the most probable configurations 

UnA**) = -^e-^ HN ^P N (da) , (1.1) 

where the normalizing factor Zn — J e~@ HN Pn, the partition function, ensures that 
(|l.ip is a probability measure, and Pn{oIu) denotes the prior distribution on Sn- 
The prior distribution is typically a product measure (see for instance Q2.2[l ) which 
describes non-interacting particle, or equivalently describes the system at infinite 
temperature ft = 0. At the ft — limit the particle interactions included in Hn 
are unimportant and thermal fluctuations, i.e., disorder, associated with the product 
structure of the prior, dominates the system. By contrast at the zero temperature 
limit, ft — > oo, interactions dominate and thermal fluctuations are unimportant; in 
this case (|1.1[> concentrates on the minimizers, also known as the "ground states", 
of the Hamiltonian Hn over all configurations a. Finite temperatures, < ft < oo, 
describe intermediate states to these two extreme regimes, including possibly phase 
transitions, i.e., regimes when as parameters, such as the temperature, change, the 
system exhibits an abrupt transition from a disordered to an ordered state and vice 
versa, or between different ordered phases. 

The objective of (equilibrium) computational statistical mechanics is the simula- 
tion of averages over Gibbs states, (jl.ip of observable quantities f(a) 

W/] = / f(^Np(da). (1.2) 

Due to the exceedingly high dimension of the integration, even for moderate values 
of the system size N, e.g., |<S/v| = 2 N for the standard Ising model, such averaged 
observables are typically calculated by Markov Chain Monte Carlo (MCMC) methods, 
|27j . Nonetheless, mesoscale morphologies, e.g., traveling waves and patterns, are 
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beyond the reach of conventional Monte Carlo methods. For this reason coarse- 
graining methods have been developed in order to speed up molecular simulations. 

We briefly discuss the mathematical formulation and numerical analysis challenges 
arising in coarse-graining of an equilibrium system described by . We rewrite the 
microscopic configuration a in terms of coarse variables r\ and corresponding fine 
variables £ so that a = (i],^). We denote the configuration space at the coarse level 
by Sm and we denote by F the coarse-graining map F : Sn — > Sm , Fcr = 77 g Sm- 
The coarse-grained system size is denoted by M, while the microscopic system size is 
N = QM, where we refer to Q as the level of coarse-graining, and Q = 1 corresponds 
to no coarse-graining. 

At the coarse-grained level one is interested in observables f(rf) which depend only 
on the coarse variable rj and a coarse-grained statistical description of the equilibrium 
properties of the system should be given by a probability measure p,M,p{drf) on Sm 
such that the average (the expected value) of such observable is same in the coarse- 
grained as well as fully resolved systems. This motivates the following definition. 

Definition 1.1. The exact coarse-grained Gibbs measure p>M,8 *s defined by 

PmA a ) = mA F ~HA)), (1.3) 

for any (measurable) set A C Sm or, equivalently, 

fWfiMAdv) = I fFWwAfo)- (1-4) 



for all (bounded) f : Sm R- 

Slightly abusing notation we will write \Xm b = I^N,B °F _1 in the sequel. In order 
to write the measure J1m,b m a more convenient form we first compute the exact 
coarse-graining of the prior distribution Pjy {da) on Sn 

P M {dv) =P N oF~ 1 , 

The conditional prior probability Pjy (da \ rf) of having a microscopic configuration a 
given a coarse configuration rj will play a crucial role in the sequel. Recall that for a 
function g(a) the conditional expectation is given by 

%N = j 9{a)P N {da\rj). (1.5) 

We now write the coarse-grained Gibbs measure P>m,b using a coarse-grained Hamil- 
tonian Hm{ti). 

Definition 1.2. The exact coarse-grained Hamiltonian Hm(v) * s given by 

e -pH M (r,) = E [ e -BH N I ^ 4 ( L g) 

This procedure is known as a renormalization group map, [181 115) . Note that the 
partition functions for Hn and Hm coincide since 



l3HN P N (da) = / / e-? HN P N (da \ r,)P M {dri) = / 6^ Sm P M {drf) = Z M • 



Hence for any function f(q) we have 

f{n) m Ada) = J f{ v )±.e^ H -P N {da) = J f(v)^ / <- ,1a j ,,)/>„ {,!,,) 



= I m^e-^^P M (dri), 
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and thus the coarse-grained measure Hm,/3 {drf) m (| 1 -3p is given by 

fiuAdv) = -^e-^Mp M (dv) • (1.7) 

Although typically Puidf]) is easy to calculate, see e.g., (|2.3|) . the exact computation 
of the coarse-grained Hamiltonian Hm(v) given by (jl.7|) is, in general, an impossible 
task even for moderately small values of N. 

In this paper we restrict our attention to lattice systems, and our main result is 
the development of a general strategy to construct explicit numerical approximations 
of the exact coarse-grained Hamiltonian Hm(j]) in the physically important case of 
combined and competing short and long range interactions. Essentially we construct 
an approximate coarse-grained energy landscape for the original complex microscopic 
lattice system in Section [2] We show that there is an expansion of Hm(t]) into a 
convergent series 

H M ( V ) = H§\rj) + H$(n) + H$(r,) + error (1.8) 

by constructing a suitable first approximation 5$ (77) and identifying small parame- 
ters to control the higher-order terms in the expansion. Truncations including a first 
few terms in (|1.8|) correspond to coarse-graining schemes of increasing accuracy. In 
order to obtain this expansion we rewrite (|1.6|) as 

H M ( V ) = H$(r)) - ~ ] gE[e-K a »- a *to» 1 77] • (1-9) 

We need to show that the logarithm can be expanded into a convergent series, uni- 
formly in N, yielding eventually an expression of the type (|1.8[) . However, two in- 
terrelated difficulties emerge immediately: (a) the stochasticity of the system in the 
finite temperature case yields the nonlinear expression in (|1.9j) which in turn will need 
to be expanded into a series; (b) the extensive nature of the microscopic system, i.e., 
typically the Hamiltonian scales as Hm = O(N), does not allow the expansion of the 
logarithm and exponential functions into the Taylor series. 

For these reasons, one of the principal mathematical tools we employ is the cluster 
expansion method, see 36^ for an overview and references. As we shall see in the course 
of this paper cluster expansions will allow us to identify uncorrelated components in 
the expected value ]SL[e~P( HN ~ H M (v)) | ^ which in turn will permit us to factorize it, 
and subsequently expand the logarithm in (|1.9[) in order to obtain the series (|1.8p . 
The coarse-graining of systems with purely long-range interactions was extensively 
studied using cluster expansions in [22] [2] [21]. Here we are broadly following and 
extending this approach. However, the presence of both short and long-range interac- 
tions presents new difficulties and requires new methods based on the ideas developed 
in |31[ [3] . Short-range interactions induce sub-grid scale correlations between coarse 
variables, and need to be explicitly included in the initial approximation ZfS (77). To 
account for these effects we introduce a multi-scale decomposition of the Gibbs state 
(|l.ip into fine and coarse variables, which in turn allows us to describe, in a explicit 
manner, the communication between scales for both short and long-range interactions. 
Furthermore, the multi-scale decomposition of ([l.ip can also allow us to reverse the 
procedure of coarse-graining in a mathematically systematic manner, i.e., reconstruct 
spatially localized "atomistic" properties, directly from coarse-grained simulations. 
We note that this issue arises extensively in the polymer science literature, |38[ 129) . 
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An important outcome of the cluster expansion analysis for the approximation of 
(11.81) is the semi-analytical splitting scheme for the coarse-graining of lattice systems 
with short and long-range interactions. Presumably similar strategies could be applied 
for off-lattice systems such as the coarse-graining of polymers. The schemes proposed 
here can be split, within a controllable approximation error, into a long and a short- 
range calculation, see (|3.27j) . The long-range part, which is computationally expensive 
for conventional Monte Carlo methods, can be cheaply simulated using the analytical 
formula given in (|3.2[) in the spirit of our previous work [35] . In this case the saving 
comes from reducing the degrees of freedom by Q = N/M and compressing the 
range of interactions. For the short-range interactions we use the semi-analytical 
formulas (I4.2[) which involve precomputing coarse-grained interactions with Monte 
Carlo simulation. However, the simulation is done for a single subdomain of three 
adjacent coarse cells. The error estimates in Theorem 13.31 also suggest an improved 
decomposition to short and long-range interactions. Indeed, they imply splitting 
and rearrangement of the overall combined short and long-range potential into a new 
short-range component that includes possible singularities originally in the long-range 
interaction, e.g., the non-smooth part in a Lennard- Jones potential, and a locally 
integrable (or smooth) long-range decaying component. 

In contrast to the splitting approach developed here that allows us to analytically 
calculate the long range effective Hamiltonian (|3.3p in Q3.27[) and in parallel carry 
out the semi-analytical step for (|4.2I) . existing methods, e.g., ([Ml [25]), employ semi- 
analytical computations involving both short, as well as costly long-range interactions. 
Thus, multi-body terms, which are believed to be important at lower temperatures, 
[14] , have to be disregarded. A notable result of our error analysis is the quantification 
of the role of multi-body terms in coarse-graining schemes, and the relative ease to 
implement them using the aforementioned splitting schemes. In Section [U we further 
quantify the regimes where such multi-body terms are necessary in the context of 
a specific example. In [2] the necessity to include multi-body terms in the effective 
coarse-grained Hamiltonian was first discussed in a numerical analysis context for 
systems with singular (at the origin) long-range interactions. 

Cluster expansions such as (|1.8[) can also be used for constructing a posteriori 
error estimates for coarse-graining problems, based on the rather elementary observa- 
tion that higher-order terms in (|3 . 33[) can be viewed as errors that depend only on the 
coarse variables n. In [20] we already employed this type of estimates for stochastic 
lattice systems with long-range interactions in order to construct adaptive coarse- 
graining schemes. These tools operated as an "on-the-fly" coarsening/refinement 
method that recovers accurately phase-diagrams. The estimates allowed us to change 
adaptively the coarse-graining level within the coarse-graining hierarchy once suitably 
large or small errors were detected, and thus to speed up the calculations of phase di- 
agrams. Adaptive simulations for molecular systems have been also recently proposed 
in |34] , although they are not based on an a posteriori error analysis perspective. Fi- 
nally, the cluster expansions necessary for the rigorous derivation and error estimates 
of the schemes developed here rely on the smallness of a suitable parameter intro- 
duced in Theorem 13.31 see Q3.32p . In Section [4] we construct an a posteriori bound 
for this quantity that can allow us to track the validity of the cluster expansion for a 
given resolution in the course of a simulation. This approach is, at an abstract level, 
similar to conditional a posteriori estimates proposed earlier in the numerical analysis 
of geometric partial differential equations, [T31 126] . 

Further challenges for systems with short and long-range interactions not dis- 
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cussed here include: error estimates for observables/quantities of interest, the de- 
velopment of coarse-grained dynamics from microscopies, phase transitions and esti- 
mation of physical parameters, such as critical temperatures. Work related to these 
directions for systems with long-range interactions have been carried out in [23], [5] 
and g]. 

The paper is organized as follows. In Section [2] we present the microscopic Ising- 
type models with short and long-range interactions and introduce the coarse-graining 
maps and the resulting coarse-grained configuration spaces. In Section [3] we discuss 
our general strategy for the analysis of systems with short and long-range interactions 
and present our main results. In Section [4] we discuss semi-analytical coarse-graining 
schemes and their applications to specific examples. Section [S] is devoted to the 
construction of the cluster expansion and to the proof of convergence of our schemes. 
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2. Microscopic lattice models and coarse-graining. We consider an Ising- 
type model on the <i-dimensional square lattice An := {x = (xi, ■ ■ ■ ,x^) e Z d ; < 
%i < n — 1} with N = n d lattice points. For simplicity we assume periodic boundary 
conditions throughout this paper although other boundary conditions can be accom- 
modated. At each lattice site x there is a spin cr(x) taking values in X = {+1, — 1}. 
A spin configuration a = {a(x)} x eA N on the lattice A at is an element of the configu- 
ration space Sn '■= S AjY . For any subset X c An we denote ax — {a{x)} x ex £ S x 
the restriction of the spin configuration to X. Similarly, for a function f : Sn — > K 
we denote fx the restriction of / to Y, x . The energy of a configuration a is given by 
the Hamiltonian 

H N (a) = H s N (a) + H l N (a) , (2.1) 

which consists of a short-range part H N and a long range part H l N . For the short- 
range part we have 

H N (a) = £ U x (a), 

XCA N 

where the short-range potential U = {Ux, X C with U x '■ ^ x — > K, is 

translation invariant (i.e., Ux+ y — Ux for all X C Z d and all y G Z d ) and has 
the finite range S (i.e., Ux = whenever diam(A) > S). We define the norm 
ll^ll = SxD{o}|diam(x)<s Halloo where the norm || • is the standard sup-norm 
on the space of continuous functions. A typical case is the nearest-neighbor Ising 
model 

H N (a) = K o-(x)*(y) , 
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where by (x, y) we denote summation over the nearest neighbors. For the long-range 
part we assume the form 

h 1 n {o-) = ~ £ E^-^W- 

where the two-body potential J has the form 

J(*-v) = j ? v(~\x-y\) , 

for some V (£ C 1 ([0, oo)). The factor 1/L d in© is a normalization which ensures that 
the strength of the potential J is essentially independent of L, i.e., J2 x ^o \ J( X )\ — 
J \ V(r)\dr. For example, if we choose V such that V(r) = for r > 1 then a spin at 
the site x interacts with its neighbors which are at most L lattice points away from x 
and in this case L is the range of the interaction J. It is convenient to think of L as 
a parameter in our model and more precise assumptions on the interactions will be 
specified later on. 

The finite-volume equilibrium states of the system are given by the canonical 
Gibbs measure and P^ida), the prior distribution on Sn, is a product measure 

P N (da) = l[ P x {da(x)). (2.2) 

A typical choice is P x (a(x) = +1) = h and P x (a(x) = —1) = i, i.e., independent 
Bernoulli random variables at each site x € A^v ■ For the sake of simplicity we consider 
Ising-type spin systems, but the techniques and ideas in this paper apply also to Potts 
and Heisenberg models or, more generally, to models where the "spin" variable takes 
values in a compact space. 

2.1. Coarse-graining. In order to coarse- grain our system we divide the lattice 
Ajy into coarse cells and define coarse variables by averaging spin values over the 
coarse cells. We partition the lattice A^ into M — m d disjoint cubic coarse cells, each 
cell containing Q = q d microscopic lattice points so that N = n d = (mq) d = M Q. The 
coarse-grained (real-space) hierarchy can be build in a anisotropic way, by replacing n, 
m, q with mult i- indexes. For example, different levels of coarse-graining in individual 
coordinate directions will be given by q = (q±,... ,qd) and the power q d would be 
interpreted as q\q2 ■ ■ ■ qd- We refrain from an unnecessary generality and assume 
that the coarse-graining is isotropic, qi = ■ ■ ■ = qd = Q- We define a coarse lattice 
Am = {k — (fci, • • • , k d ) € 7L d ; < fej < m — 1} and we set A N = U kE * C/. where 
Cfe = { x € Am ; fcjg < Xi < (fej + l)q}. Whenever convenient we will identify the 
coarse cell Ck in the microscopic lattice A^r with the point k of the coarse lattice 
Am- For any configuration a k = ac k on the coarse cell Ck we assign a new spin value 

V(k) = £ o-(x) 

which takes values in S = {— Q,— Q + 2, . . . ,Q}. We denote the configuration space 
at the coarse level by Sm = S Am and we denote by F the coarse-graining map 



F : S N -> S M , cr = {a(x)} xeAN (->• r) = {v(k)} ke K 
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which assigns a configuration 7/ on the coarse lattice Am given a configuration a on 
the microscopic lattice Ajy- 

The exact coarse-grained Gibbs measure is defined in (|1.3[) for arbitrary Gibbs 
states having the form (11.7[) . Since rj(k) depends only on the spins cr(x), with x G Ck, 
the coarse-grained measure Pm is a product measure 

P M {dr,) = P N oF- 1 = H P k {dr,{k)). (2.3) 



For example if P x is a Bernoulli distribution then Pk(rj(k)) = (j}(A;)+q) (|) • Simi- 

2 

larly, we define the conditional probability measure P/v {da | 77) of having a microscopic 
configuration a on A^v given a coarse configuration r\ on Am- This measure plays a 
crucial role in the sequel since it factorizes over the coarse cells 

P N {da\rj) = H P k (da k \ v (k)), (2.4) 
fceA M 

where P k {da k \ r){k)) is the conditional probability of a microscopic configuration a k 
on Ck given a coarse configuration rj{k). 

3. Approximation strategies for Hm(ii)- I n this section we present a gen- 
eral strategy for constructing approximations of the exact coarse-grained Hamiltonian 
Hm{v) i n (|f -71) . We show how to expand HmIii) into a convergent series (If .8[) by 
choosing a suitable first approximation _ff M ^(ry) and identifying small parameters to 
control the higher-order terms in the expansions. The basic idea is to use the first 
approximation .ffjS (77) in order to rewrite (If .6[) as (| 1 .9[) . We show that the logarithm 
can be expanded into a convergent series, uniformly in N, using suitable cluster ex- 
pansion techniques. We discuss in detail the case d = 1 in order to illustrate general 
ideas in the case where calculations and formulas are relatively simple. The general 
d-dimensional case is discussed in detail in Section [5] 

We recall that the Hamiltonian Hm{<j) — H l N (a) + Hfj(a) consists of a short- 
range part H^(a) with the range S and a long-range part H l N (a) whose range is L. 
We choose the coarse-graining level q such that 

S < q < L. 

There are two small parameters associated with the range of the interactions 

S q 
e s oc — , and ei oc — . 
q L 

The first approximation is of the form 

H^=HT+H^\ (3.1) 

and two distinct separate procedures are used to define the short-range coarse-grained 
approximation Hl'j °\ as well as its long-range counterpart H l ^°\ Due to the non- 
linear nature of the map induced by (|1.9[) it is not obvious that (|3.1[) will be a valid 
approximation, except possibly at high temperatures, when /3 « 1. This fact will 
be established for a wide range of parameters in the error analysis of Theorem 13.31 
and in the discussion in Section 01 provided a suitable choice is made for 3^ an d 

rr*,(0) 
n M ■ 
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3.1. Coarse-graining of the long-range interactions. We briefly recall the 
coarse-graining strategy of [55J for the long-range interactions. Since the range of the 
interaction, L, is larger than the range of coarse-graining Q a natural first approxi- 
mation for the long-range part is to average the interaction J(x — y) over coarse cells. 
Thus we define 

= E[H l N \rj\, (3.2) 

and an easy computation gives 

#m%) = ~E £ J (MM*M0-5 E J(k,k)( V (k) 2 -Q), (3.3) 
fceAM l ^= k keA M 

where 

A simple error estimate (see [55] 12] for details in various cases) gives 



H l N (a) = H l ^ 0) (F(a)) + e L with e L = jV0(||| WIU) . 



Using this definition of H l A ^ we obtain 



-^WiV(da I r?) = e-^We^ ^W-^'')] P N (da \ r,) , (3.4) 



n (1+4)^(^1??), (3.5) 



where 



Due to the fact that P/v (da 1 77) has a product structure one can rewrite (13.51) as a 
cluster expansion, [22] (see also Section[SJ), as in (jl.8[> . The key element in that cluster 
expansion is the "smallness" of the quantity 

\J(x-y)- J(fc,0|<2-|- sup |W(z'-j/)l, (3.7) 

^ x'ec,., 



which yields asymptotics 



/jfc-^^HV^loo). (3.8) 



The estimate (|3.7p follows from regularity assumptions on V and the Taylor expansion. 

3.2. Coarse-graining of short-range interactions. For the short-range part, 
using that S < q, we write the Hamiltonian as 

H' N (a)= E H k(°)+ E W k , k+1 (a), (3.9) 
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where 

#fe» = E ^W, W k , k+l {a) = Y, Ux ^> 

XdC k Xr\C k ^$ ,XnC k+1 ^$ 

i.e., Hf. is the energy for the cell C k which does not interact with other cells, i.e., 
under the free boundary conditions, and W k . k +i is the interaction energy between the 
cells C k and C k +i- Note the elementary bound 

supW M+1 (ff)~£3 ,, - 1 ||Z7||. (3.10) 
(j 

The most naive coarse-graining, besides of course developing a mean-field-type ap- 
proximation, consists in regarding the boundary terms W k , k +i as a perturbation. We 
have then, formally, 

= e -m' M m ( V )+e L +e s TT e -PU s k m ( Vk ) ) 

fee Am 

where the one-body potential 

^' (0) fe) = -ilog | e-^^P fe (da fe h(fc)) 

is the exact coarse-grained Hamiltonian for the cell C k with free boundary conditions. 
As a result an initial guess for the zero order approximation could be 

HT(v)+Y,U S k M (Vk). (3.H) 

k 

However, this approach appears to be rather simplistic in general since the correlations 
between the cells induced by the short-range potential have been completely ignored. 
While this approximation may be reasonable at high temperatures it is not a good 
starting point for a series expansion of the Hamiltonian using a cluster expansion. 
Instead we need to adopt a more systematic approach outlined in the next section. 

3.3. Multiscale decomposition of Gibbs states. This approach provides the 
common underlying structure of all coarse-graining schemes at equilibrium including 
lattice and off-lattice models. It is essentially a decomposition of the Gibbs state Ijl.ljl 
into product measures among different scales selected with suitable properties. We 
outline it for the case of short-range interactions where we rewrite the Gibbs measure 

CEZED as 

mAd°) ~ e-P H »^P N (da) = e- pH "^P N {da \ r,)Pu(dv) ■ 

We use the notation ~ meaning up to a normalization constant, i.e., in the equation 
above we do not spell out the presence of the constant Zn ■ We now seek the following 
decomposition of the short-range interactions 

e- pa '"^P N (da | v) = R{v)A(a)v(da\ri) , (3.12) 



where 
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(a) R(ri) depends only on the coarse variable r\ and is related to the first coarse-grained 
approximation H^' ^ via the formula 

R( v ) = e -^M (0) W , A{a)v{da\ri) = e^^"^^ (7l) ) P N {da \ rj) , (3.13) 

(b) A(a) has a form amenable to a cluster expansion, i.e., for d = 1 

for some K C A^/. The function is small and moreover ^-(cr) depends on the 
configuration c only locally, up to a fixed finite distance from Ck- In the example at 
hand (for d = 1) we have $fc(cr) = $fe(cr fc_:L , cr fc+1 ). 

(c) The measure v(da\rf) has the general form 

v{da\i)) = ]J i/ fc (db-|f|) , (3.15) 

fee Am 

where v k (da \ rf) depends on a and r\ only locally up to a fixed finite distance from 
Cu- In the example at hand Vk{d<J \ rf) depends only on the configuration on Ck-i U 
Ck U Ck+i- Even though the measure v(da\r]) is not a product measure, the fact that 
this measure has finite spatial correlation makes it adequate for a cluster expansion, 
see (|3.26p and Section [5] 

Although here we described the multiscale decomposition of the Gibbs measure 
for the case of short-range interactions, the results on the long-range interactions, 
discussed earlier, can be reformulated in a similar way. In particular, Q3.4p and (|3.4p 
can be rewritten as 

e -!3H' N (*)p N ( d(J | ^ = R ( K r ) )A(a)v(da\r 1 ) , (3.16) 
where R(rj) = e~^M W, v(d<r\r\) = Pjy(da \ rj), and 

A(*)=e^( H '^-^r(v)) = Y[ (1 + 4). (3.17) 

j,keA M 

We recall that in analogy to (|3.15[) , the product structure of v(da\ri) — P N (da\r]) 
allows us to carry out a cluster expansion for the long-range case, and obtain a con- 
vergent series such as (|1.8p . thus yielding an expansion of the exact coarse-grained 
Hamiltonian H l M , [22] . 

We note that (|3~T2|) . used here as a numerical and multiscale analysis tool in order 
to derive suitable approximation schemes for the coarse-grained Hamiltonian, was first 
introduced in [301 EU [3] for the purpose of deriving cluster expansions for lattice sys- 
tems with short-range interactions away from the well-understood high temperature 
regime. 

3.4. Coarse-graining schemes in one spatial dimension. We sketch how 
to obtain a decomposition such as (|3.12|) for d = 1 and construct suitable R(rj). We 
split the one-dimensional lattice into non-communicating components, for instance, 
even- and odd-indexed cells and write 

e- pB »P N (da\ri) = J[ \e- 0i - Wk -^ +Wk - k ^e- f3HS >'P k {dcT k \ n{k))] x 

k: odd 

[] e-^P k (da k \r](k)). (3.18) 

k: oven 
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In (|3.18[) we will normalize the factors for k odd by dividing each factor with the 
suitably defined corresponding partition functions for the regions C k and Cfc_i UC k U 
Ck+i- 

Definition 3.1. We define the partition function with boundary conditions <r k ~ 1 
and o~ k+1 , i.e., 

Z k (n{k)-a k -\a k+1 ) = [ e-^ Wh -^ k+Wk ' k +^e-^P k {da k | . (3.19) 



In order to decouple even and odd cells we define the partition function with free 
boundary conditions on C k -\ and boundary condition a k+1 on Cfe+i, i.e., 

Z k {n{k)-0,a k+1 ) = J e-^ w ^e-^P k (da k \n(k)), (3.20) 

and similarly Z k (n(k); a k ~ 1 , 0), as the partition function with free boundary conditions 
on Cfc+i and boundary condition a k ~ 1 on C k -i- We also denote by Z k (r/(k);0,0) 
the partition function for C k with free boundary conditions. We define the three-cell 
partition function with free boundary conditions 

Zfc_i 1 fc lfc+ i(»7(fe-l),»?(fc),»7(fc + l);0,0) = 

e -0(#2_ 1 +w fc ,*_i+irj+w fc ,* + i+fl2 +1 ) x 
Pk-^da*- 1 I n(k - l))P k (da k | n(k))P k+1 (da k+1 \ V (k + 1)) . (3.21) 

The key to the decomposition and eventually to the cluster expansion is the intro- 
duction of a "small term" analogous to (|3.8p . 

Definition 3.2. 

h-x, k+ im),o- ,a ) - z^ky^^z^ik)-,^^)- 1 (3 - 22) 

An important element in the cluster expansion in Section [3] is the estimation of the 
terms k+1 . However, a straightforward estimate based on (|3.10p would yield 

fLi,k+Mk);a k -\a k+1 ) ~ /3S\\U\\ . (3.23) 

We rewrite 



Z k (ri(ky,a k -\<T k + 1 )=(f k ^ tk+1 (T,(ky,a k - 1 ,a k+1 ) + l) x 

Z k (T t (k);0,<T k + 1 )Z k (r](k);a k - 1 ,0) 
Z k (ri(k); 0,0) 



(3.24) 



In (|3.18p we now divide and multiply each factor with k odd by Z k (a k 1 ,cr +1 ) 
and use the formula (13.24)) . Furthermore, we multiply each factor with even k by 
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£fc-i,k,k+i(0,0) and obtain 
e -P H »P N (da\r)) = 

n ^(o,o)- x n ^-i,k,k+i(o,o) n + x ) x ( 3 - 25 ) 

fc: odd k: even fe: odd 



flfa) = A(<r) 



n e v i a k i n \\ TT e fcZ fc+n cr )UJZ fc _i(U,tr J . fe . 

^-spw-'W' l,W) fc n — — p * 

= K^M) 

(3.26) 

where we have used that 

Z k (0,a k+1 )Z k (a k -\0)= JJ Z fc+ i((7 fc , 0)Z fc _i(0, a fe ) . 

fc: odd k: even 

It is easy to verify that v(da \ rf) defined in (|3.26[) is a normalized measure and has the 
form required in condition (c) of the multiscale decomposition of the Gibbs measure. 
The factor R(rj) defined in (|3. 251) gives the first order corrections induced by the 
correlations between adjacent cells. Putting together the analysis for short and long- 
range interactions we obtain the main result formulated as a theorem. 
Theorem 3.3. Let 

H^(V) = HT(V) + Hlf\v) (3-27) 
where H 1 ^ ■ (rj) is given in (|3.2[) and (I3.3[) and 

#m%) = E U s k M (v(k))+ E + (3.28) 

k: odd k: even 

with the one-body interactions 

U S k i0) (v(k)) = ~ log Z k ( V (k); 0,0), (3.29) 
and the three-body interactions 

^-! w (#-i)#).# + i)) = 

~ log Z k ^ k+1 ( v (k - 1), V (k), r,{k + 1); 0, 0) , (3.30) 

where Z k and ^k-i,fe,fe+i are given by US. 20\) and \3.21\i respectively. Then 
1. we have the error bound 

\ q L 

for a short-range potential with the range S « q « L. The loss of information when 
coarse- graining at the level q is quantified by the specific relative entropy error 

| w „F-.) = (ffl + i&). (3.3!) 
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2. There exist S a > and Si > such that if 

sup sup |/fc_i,fc+ifa(fc);0- fc-1 ,CT fc+1 )| < S Q , sup sup \f jk (a j ,a k )\ < Si , 

(3.32) 



(0) 
M 



where fc+1 and /L are <?u>en &?/ VS. 22^ and VJ. 61) respectively, then Hm — H 
expanded in a convergent series in the parameter 8 ~ (MZM + 9/3]|W]Uc ^ 

ff*f(l) - ff£ } (»7) + ff&'fa) + ' ■ ' + ^ } (»7) + MO(SP +1 ) . (3.33) 



Remark 3.1. The error estimate (|3.3ip suggests qualitatively an estimate on 
the regimes of validity of the method, and on the "optimal" level, q = q pt, when we 
restrict to the regime S < q < L, where S and L are the respective interaction ranges 
for short and long-range potentials. The corresponding error is then 




(3.34) 



The application of Theorem [3?3] requires to check the validity of (|3 . 32|) . Certainly 
the conditions (13.81) and (|3.23[) are satisfied in suitable regimes, see also Section [5] for 
more details. More interestingly, for specific examples these conditions can be verified 
directly, we refer to Section 2J In particular, in (|4.9p and (14. 13)) we even obtain an 
upper bound that depends only on the coarse observables. This allows us to check the 
conditions (|3.32[) (dictated by the cluster expansions) computationally in the process 
of a Monte Carlo simulation involving only the coarse variables 77. 

On the other hand, in [30j [31] , the short-range condition in (|3.32[) is taken as 
an assumption. In one dimension, this condition holds up to very low temperatures 
while in dimension d > 2 this condition can be satisfied in the high-temperature 
regime, see for example the analysis in [3] where similar conditions are used for the 
nearest-neighbor Ising model in the dimension d = 2 all the way up to the critical 
temperature. 

Finally, we note that a similar strategy to coarse-grained short and long-range 
interactions can be used in any dimension, as we discuss in Section [5] In the multi- 
dimensional case we split the domain into boxes of size larger than the range of the 
interaction so that the next-to-nearest coarse cells are independent. In one dimension, 
this procedure gives rise to the separation into odd- and even- indexed coarse cells, 
while in higher dimensions it is done in a recursive manner, proceeding one dimension 
at a time. Then by freezing the configurations on the collection of independent coarse 
cells (resulting to the one-body coarse-grained terms) we create further correlations 
which couple the remaining cells. This fact in one-space dimension yields the three- 
body terms, noting that possible two-body coarse-grained correlations are contained 
therein, see also (|4.8p . We also remark that coarse-graining schemes for the nearest- 
neighbor Ising model, involving only two-body interactions were recently proposed in 

Outline of the proof: Using the coarse-grained approximation H$ (77) the decompo- 
sition p.!2|) can be rewritten as R(rj) = eT* 1 ** and thus we obtain 

H M (rj) = H$(r,) - i log J A(a)v{da\r,) , 
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where A, and v„ are given abstractly in A3. 13)) and are denned both for short and 
long-range interactions in analogy to (I3.26|) . The construction of the series in (I3.33[) 
relies on the cluster expansion of the type 

A(a) ee Ha + 4.) n (1 + //_ M+1 ) = n /« ( 3 - 35 ) 

i<j i-odd Geg M {i,j}eE(G) 

where 

fh or fjj , if i even and j = i ± 2 

/■j otherwise, 

and 5m is the set of all graphs on M vertices, where M is the total number of coarse 
cells. Such an equality and the complete proof is carried out in Section [5j In turn, 
the terms on the right hand side of (|3.35j) give rise to the expansion (|3.33j) and the 
corresponding higher-order corrections. 

3.5. A posteriori error estimates. In (52] we introduced the use of cluster 
expansions as a tool for constructing a posteriori error estimates for coarse-graining 
problems, based on the rather simple observation that higher-order terms in p.33[) 
can be viewed as errors that depend only on the coarse variables rj. Following the 
same approach an a posteriori estimate immediately follows from (|3.33p . 

Corollary 3.4. We have 

I W,fi o F- 1 ) = m m [S(r,)] + log (E ( o, [e-^]) + 0(S 2 ) , 

where the residuum operator is S(rj) — i?U (77). 

In [20] we already employed this type of estimates for stochastic lattice systems 
with long-range interactions, in order to construct adaptive coarse-graining schemes. 
These tools operated as an "on-the-fiy" coarsening/rehnement method that recov- 
ers accurately phase-diagrams. The estimates allowed us to change adaptively the 
coarse-graining level within the coarse-graining hierarchy once sufficiently large or 
small errors were detected, thus speeding up the calculations of phase diagrams. Ear- 
lier work that uses only an upper bound and not the asymptotically sharp cluster 
expansion-based estimate can be found in [5] [7] . 

3.6. Microscopic reconstruction. The reverse procedure of coarse-graining, 
i.e. reproducing "atomistic" properties, directly from coarse-grained simulation meth- 
ods is an issue that arises extensively in the polymer science literature, [3H [53] . The 
principal idea is that computationally inexpensive coarse-graining algorithms will re- 
produce large scale structures and subsequently microscopic information will be added 
through microscopic reconstruction, for example the calculation of diffusion of pene- 
trants through polymer melts, reconstructed from CG simulation, [29] . 

In this direction, the CGMC methodology discussed in this section can provide a 
framework to mathematically formulate microscopic reconstruction and study related 
numerical and computational issues. Indeed, the conditional measure A(a)u(da\r)) in 
the multi-scale decompositions (|3.12l) and (j3. 16[) can be also viewed as a microscopic 
reconstruction of the Gibbs state (|1.1[) once the coarse variables ij are specified. The 
product structure in (|3.14p and (|3.15[) allows for easy generation of the fine scale 
details by first reconstructing over a family of domains given only the coarse-grained 
data and gradually moving to the next family of domains given now both the coarse- 
grained data and the previously reconstructed microscopic values. 



fij — 
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In view of of this abstract procedure based on multiscale decompositions such as 
(|3.12p . we readily see that the particular product structure of the explicit formulas 
(|3.25[) and (|3.26p for the case of the dimension d = 1 yields a hierarchy of recon- 
struction schemes. A first order approximation can be based on the approximation 
A(a) ~ 1 (cf. (13331) . (prig) ): 

(a) first, R(r]) defined in Q3.25P provides the coarse-graining scheme, which will 
produce coarse variable data rj(k) for all k; 

(b) next, we reconstruct the microscopic configuration er cvcn consisting of the 
er fc 's in all boxes (coarse-cells) with k even using the measure i/k(d<r\i]) := 

^ z^i^fc+i^o) 1 ^ ^ " ^ >k ^ (jk I^W)' conditioned on the coarse configura- 
tion rj{k) from (a) above; 

(c) finally, we reconstruct the microscopic configuration in the remaining boxes 

~p( H k + W k-l k + W k k + l) 

with k odd using i/ k (da\i]) := $ g^^r^i+ry Pk(do k \ f](k)), given the 

coarse variable rj(k) from step (a), and the microscopic configurations <j even 
from step (b). 

We note that this procedure is local in the sense that the reconstruction can be carried 
out in only the "subdomain of interest" of the entire microscopic lattice A^r; this is 
clearly computationally advantageous because microscopic kMC solvers are used only 
in the specific part of the computational domain, while inexpensive CGMC solvers 
are used in the entire coarse lattice Am- 

Further discussion on the numerical analysis issues related to microscopic recon- 
struction for lattice systems with long-range interactions can be found in [21] [37l [23l 
[24]. 

4. Semi-analytical coarse-graining schemes and examples. Next we dis- 
cuss the numerical implementation of the effective coarse-grained Hamiltonians de- 
rived in Theorem 13.31 We begin with a general implementation scheme and we sub- 
sequently investigate further simplifications for particular examples in one space di- 
mension. 



4.1. Semi-analytical splitting schemes and inverse Monte Carlo meth- 
ods. One of the main points of our method is encapsulated in (|3 . 2T|) : the compu- 
tationally expensive long-range part for conventional Monte Carlo methods can be 
computed by calculating the analytical formula given in (|3.2[) in the spirit of our pre- 
vious work [22j . Then we can turn our attention to the short-range interactions where 
Monte Carlo methods, at least for reasonably sized domains, are inexpensive. More 
specifically for the evaluation of the short-range contribution in (|3.27l) we introduce 
the normalized measure 

where the sum is computed with free boundary conditions on C k and Zk(r)(k); 0, 0) is 
accordingly defined as in (|3.20l) . Thus (|3.28[) can be rewritten as 

H s A f =Y,U s k m (v(k))+ £ V£l C ? M+1 fa(*-l), V (k),r,(k + 1)), (4.2) 

fcGA fe: even 
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where, based on (|3.28l) and (|4.ip , we denned the three-body coarse interaction poten- 
tial 

xft-ifd^ 1 | v(k - l))Pk(da k \ V (k))P k+1 (da k+1 \ V (k + 1)) . (4.3) 

The main difficulty in the calculation of (|4.3p is that for the three-body integral one 
needs to perform the integration for all possible combinations of the multi-canonical 
constraint. On the other hand all simulations involve only short-range interactions 
and need to be carried out only on three coarse cells, rather than the entire lattice. 
Practically, the calculation of (|4.3j) can be implemented using the so-called inverse 
Monte Carlo method, [25] . We sample the measure Pk using Metropolis spin flips 
and subsequently we create a histogram for all possible values of rj{k) — Yl x eC <r(x). 
Then we compute the above integral by using the samples which correspond to the 
prescribed values rj(k — 1), rj(k) and rj(k + 1). 

A complementary approach in order to further increase the computational effi- 
ciency of the schemes presented in Theorem 13.31 is to rearrange the splitting based 
on the size of the error in (|3.31|) . Indeed, these estimates suggest a natural way to 
decompose the overall interaction potential into: (a) a short-range piece J s including 
possible singularities originally in J, e.g., the non-smooth part in the Lennard- Jones 
potential, and (b) a locally integrable (or smooth) long-range decaying component, 
J;. Thus, if K(x,y) is the short-range potential in (|2.1j) we can rewrite the overall 
potential as 

K(x, y) + J(x, y) = J s (x, y) + Jj(x, y) . (4.4) 

In this way the accuracy can be enhanced by implementing the analytical coarse- 
graining (|3.3|) for the smooth long-range piece Ji(x, y), and the semi- analytical scheme 
(I3.28P for the "effective" short-range piece J s (x,y). 

Remark 4.1. Existing methods, e.g., [2], employ an inverse Monte Carlo com- 
putation involving both short and long-range interactions, and due to computational 
limitations have to disregard multi-body terms such as the ones considered in the 
method proposed here. The splitting approach developed here allows us to calculate 
analytically the approximate effective Hamiltonian for the costly long-range interac- 
tions, (|3.3|) in (|3.27l) or (I4.4[) . and in parallel carry out the inverse Monte Carlo step 
for (|4.2I) . The necessity to include multi-body terms in the effective Hamiltonian was 
first discussed in [2 together with their role in the proper coarse-graining of singu- 
lar short-range interactions. We further quantify the regimes where such multi-body 
terms are necessary in the context of a specific example. 

4.2. A typical example: improved schemes and a posteriori estimation. 

We examine the derived coarse-graining schemes in the context of a specific, but rather 
typical example. We consider the Hamiltonian 

H N {a) = H s N {o-) + H l N (a) := K ]T a(x)a(y) - \ £ J(x - y)a(x)a(y) (4.5) 

(*,y) (x,y) 

where by (x,y) we denote summation over the nearest neighbors, i.e., \x — y\ = 1, and 
by (x,y) the long range summation as in @. Although we follow the splitting strat- 
egy discussed in the previous paragraph we present a simplified numerical algorithm 
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by carrying out further analytical calculations. Not surprisingly, such calculations 
allow not only for easier sampling in the semi-analytical calculations of the inverse 
Monte Carlo, but give additional insight on the nature of multi-body, coarse-grained 
interactions. 

For the short-range contributions, given a coarse cell C k with q lattice points, we 
denote by x\, . . . ,x q the lattice sites in C k . With this notation, following (|4.3p the 
short-range three-body interaction is given by 

xA-if^" 1 I V(k - l))P k (da k | V (k))P k+1 (da k+1 \ V (k + 1)) . (4.6) 

The main difficulty in computing the second term is the conditioning on the 
coarse-grained values r](k — l),T](k),r)(k + 1) over three coarse cells. At first glance 
this requires to run multi-constrained Monte Carlo dynamics for every given value of 
the 77's, i.e., for q 3 variables. However, as we show in the sequel, when dealing with 
a particular example, e.g., the nearest neighbor interactions, the computationally 
expensive three-body term reduces to product of one-body terms. We first rewrite 

where we set 

a = cosh(/3K) , b = smh((3K) , A = taah(PK) . 
Moreover, we introduce the one- and two-point correlation functions 

*k(%):= / a{x)P k {da k \r](k)) and (%):=/ a(x)a{y)P k (da k \ v (k)) . 



By symmetry we have that Q^ 1 = and similarly, consider 9 for x = xi and 
y = x q . Furthermore, these functions depend on k only via the coarse variable rjk, 
hence we now define 



$ x ( % ) := / a(xt)P k (da k \ r,(k)) and <f> 2 (r] k ) := / a( Xl )a(x q )P k (da k \ rj(k)) . 

(4.7) 

It is a straightforward computation to show that 

V k s %, k+ Mk ~ lMk),ri(k + 1)) = -~ log a - 

-ilog (l - A$ 1 (r ? (fc - l))$ 1 ( ?? (fc)) - A$ 1 (r ? (fc))$ 1 (77(fc + 1)) 

+A 2 $ 1 (r/(fc - l))$ 2 (?7(fc))$ 1 (??(fc + 1))) (4.8) 

Although these are three-body interactions, the additional analytical calculations re- 
duce their computation to the nearest-neighbor Monte Carlo sub-grid sampling of 
(|4.7[) . Moreover, from f|3 . 22[) we have 



fS , m k-i k + i, \ 2 a k ' 1 (x q )a k+1 (xi ) [$ 2 (rj(k) 

/fc-i,fc+iW ho , ) ^ 1 _Xak-i( Xq )^(r](k)))(l-Xa k + l {x 1 )^(r](k))) 
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thus the following estimate holds for some C > 

sup \fu,k+i\ ^ CA2 l $2 W fc )) ~ l^(v(k))?\ = e(w A) , (4.9) 

where the right-hand side G is an a posteriori functional in the sense that it can be 
computed from the coarse-grained data. In fact, we can estimate the a posteriori error 
indicator by an analytical formula. A high temperature expansion yields 

$!(%) = E[<r(aO 1 n] + 0(A) = ^ + 0(A) (4.10) 

q 

<£ 2 (r, k ) = E[a(x)a(y) | rj\ + 0(A) = + 0(A) . (4.11) 



Then. 



Thus the validity of Theorem 13.31 and the derived coarse-grained approximations can 
be conditionally checked during simulation by 

2 / „2 



0( Vk ; A) ~ A 2 |$ 2 - ($i) 2 | = A 2 g f (g _^ } + 0(A 3 ) . (4.12) 



A / 7] 



sup \!Li, k+ i\<C )+0{\ A ). (4.13) 

We note that (|4.13p suggests a quantitative understanding of the dependence of 
the coarse-graining error for the nearest-neighbor Ising model. The error increases, 
(a) when the parameter A 2 increases, i.e., at lower temperatures/stronger short-range 
interactions, (b) when the level of coarse-graining q decreases, and (c) at regimes where 
the local coverage r\ is not uniformly homogeneous, i.e., away from the regime ry » ±g. 
Such situation occurs, for example, around an interface in the phase transition regime. 
This is the case even in one dimension if long-range interactions are present in the 
system. 

5. Proofs. In this section we first construct and prove the convergence of the 
cluster expansion. We formulate the proofs in the full generality assuming a d- 
dimensional lattice. Thus coordinates of lattice points are understood as multi-indices 
in Z d . We start by constructing the a priori coarse-grained measure induced by 
the short-range interaction. We perform a block decimation procedure following the 
strategy in [21] and partition Am into 2 d -many sublattices of spacing 2q. Let e a , 
a = 2, 3, . . . , 2 d be vectors (of length q) along the edges of Am as demonstrated in 
Figure E3] for d = 3. We write the coarse lattice as union of sub-lattices 

Am - U 2d =1 A M , (5.1) 

where K\j = 2 Am, A m = A M + e 2 and A^ 1 = A% f + e a+ i, for a = 1, . . . , 2 d - 1. 
Given a coarse cell Ck we define the set of neighboring cells by 

8Ck ■■= U{; : |||_fc||=l}Cj| , 

where \\l — k\\ := max i=lr .. jd \l t — fc, ; |. We also let D k :— C k U dC k - 

Given a sublattice A^ we denote by a a the microscopic configuration in all the 
cells Ck G and by a >a the configuration in K^ M for all /3 > a. We also define a 
function p : Am — > {1, • ■ • , 2 d } such that for fc £ Am, we have p(k) = a if Ck £ A M . 
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Fig. 5.1. The sublattices A M covering the coarse lattice Am- The vectors e a defining transla- 
tions of the first sublattice A M are depicted for d = 2,3. The cells on the two-dimensional lattice 
are numbered with values of a = 1, . . . , 4 according to what sublattice A?» they belong. 



We split the short-range part of (|2.1[) 

*W=£ E E W k {o»;v> a ), 

where, for fc £ A^ f , the terms Hk{cr a ) are the self energy on the boxes Ck given by 

m°°)= E mo- 

xcc fc 

Moreover, the energy due to the interaction of Ck with the neighboring cells is given 

by 



W k (a a ;a >a )= £ U x (a a V a >a ) 



XCD k 

where a a V<r >Q is the concatenation on A^ and A^ Q . Now we construct the reference 
conditional measure v n under the constraint of a fixed averaged value rj — {i]{k)}keA M 
on the coarse cells. 

Step 1. The starting point is a product measure on Ck for k g Aj^. We let A\{k) = Ck 
and after appropriate normalization we obtain 

IJ ^(fc);^ 1 ;^)) ^(da 1 ) 



(5.2) 



where 



k£A\ 



Z(Ai(fc);CT >1 ;?7) 



(5.3) 



is the new prior measure on with boundary conditions a >x and the canonical 
constraint i](k), k G Aj^. The partition function 

Z{A 1 (k)-a >1 -T 1 {k))= f e-x^e-W^'^Pkida 1 ) 
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depending on the boundary conditions a >1 on the set dAi(k) couples the configura- 
tions in Ci with I £ dA 1 (k). In particular, it couples the configurations a 2 and gives 
rise to a new interaction between them for which it will be shown that it is small due 
to the condition 15.11 

Step 2. Moving along the vector e 2 we seek the measure Vy 2 on {+1, —l} UkeA M Ck . 
Given the partition function Z(Ai(k); cr >1 ; i](k)) we denote by S ke2 Z the partition 
function on the same domain A\{k) as but with new boundary conditions which 
are the same as Z in the +e2 direction, free in the — e 2 and unchanged in all the other 
directions. Similarly, we denote by S ke2 Z the partition function with free boundary 
conditions in the direction +e 2 and by S® e2 Z with free boundary conditions in both 
±e 2 directions. With these definitions we have the identity 

Z{A 1 {k) ] a>\r,{k)) = iS ^ )iS ^ Z \ l + H), (5.4) 

where we have introduced the function $J. which contains the interaction between the 
variables cr >1 , and it is given by 

= Z(A 1 (k);a> 1 ;r,(k))(S kiea Z) i 



In this way we split the partition function Z into a part where the interaction between 
the cells C'k~ e2 and Ck+e 2 lS decoupled and an error part which is to be small. The 
terms in the second product contain all possible interactions in the set 

A 2 (k) = GVe 2 U C fe U C k+e2 (5.5) 

for k £ h? M with the corresponding partition function being given by 

Z{A 2 {k);a> 2 ;r 1 {k)) = J e-^^) e ^ fc ^^ >2 )(5+_ e2 ^Z)(^ +e2!e2 Z)P fe (d ( T 2 ) 

all due to the condition 15.11 



The next step is to index the new partition functions {S ke2 Z) and (S k e2 Z) (which 
are functions of a 2 indexed by k £ Aj^-) with respect to k £ h. 2 M . We have 

n (st, 2 zm e2 z)= n (ste 2 , e ^m +e2>e2 z). 

feeAJ, keK 2 M 

Then if we neglect for a moment the error term (1 + in order to define v 2 >2 we 
have to deal with the following terms 

fceA^ keA 2 M 

The terms in the second product contain all possible interactions in the set A 2 (fc), 
given in (|5.5[) for k £ K 2 M with the corresponding partition function being given by 

Z(A 2 (fc);a> 2 ;?? (fc)) = / e -^^)e-^^^ >2 )(5+_ e2ie2 Z)5 fc+e2!e2 ^)P fc (d<7 2 ). 
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By normalizing with this function we obtain the measure 

1 



v 2 >2 (da 2 )= J] 



ke\{ 



Z(A 2 (k);a> 2 ; V (k)) 



(5.6) 



Note that the factor (S^. e2 Z) _1 depends on r\ as well as on a >2 and hence we will 
need to further split it when we define a measure on the variables on which it depends. 
Summarizing the first two steps we have obtained that the left hand side of (|5.2jl is 
equal to 



[] Z(A 2 (k);a> 2 ; V (k)) U ffi^Z)- 1 J] i 1 + *l) 



keAl 



If we are interested in the case d = 1, this would be the final expression. However, 
for higher dimensions we need to repeat the above steps. We give one more step 
in order to obtain more intuition on the relevant terms and then we give the final 
expression in agreement with the result in |31) . The proof of the general formula is 
done with a recurrence argument on the number of steps and for the details we refer 
to 

Step 3. To proceed in the next step along direction e 3 we split Z(A 2 (k); <r >2 ; i](k)) 
(which couples the configurations in Ck with p{k) = 3) in the same fashion as before. 
We have 



Z(A 2 (k);a> 2 ; V (k)) = 



where 



Z(A 2 (k);a> 2 ; V (k))(Sl e3 Z) 



We further change the indices in such a way that they are expressed with respect to 
k G Aj^ and then we glue the partition functions on A 2 {k — 63) and A 2 (k + 63). 
We define 

A 3 (k) := C k L)A 2 {k- e a ) U A 2 {k + e 3 ) , 



and 

Z(A 3 (k);a >3 ; V (k)) 
The corresponding measure is 

4 3 (^ 3 )= n 



keA 3 



-Ht„-W k (a 3 -a> 3 ; V ) 



1 



(Ste 3 ,e 3 Z)(S^ e3 Z)P k (da 3 ). 



Z(A 3 (k);a>3 ]V (k)) 



k — e 3 ,e; 



z)(s, 



fc+e3 ,e3 



Z)Pk(da 3 ) 



(5.7) 
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and the left hand side of (|5.2j) is now equal to 

n [e-^e-^^p^)] n z(Mk);* >3 -Mk)) n KesZ)- 1 * 

n Ke^r 1 n n a+^^^^Ki^ 1 )- 

feeA^ keKl, keAl 

As in the previous steps we need to perform the usual actions on the partition function 
Z(A 3 (k); er >3 ; rj(k)) which will give rise to a new element A 4 (fc) with k £ A M and new 
error terms <fr£ with /c ^ A M . Furthermore, a similar splitting has also to occur for the 
factor (S® e2 Z)^ 1 which also depends on cr 4 , since the zero boundary condition involves 
only the direction e^- Related calculations will involve all the terms of similar origin 
as long as we move to new sublattices A M , with a > 4, depending on the dimension. 
Example: 2D lattice The leading term in the approximation of the coarse-grained 
Hamiltonian consists of terms that refer to four different types of multi-cell in- 
teractions 



HT= E ^gZ(A 4 (k))- E ^Z(A 4 (k)) 



r '"" " E 

k£A\, keA 2 , 



+ J2 bgZ(Mk))- E logZ(A A (k)), 

keAl, keAl, 

where A±(k) is a collection of coarse cells centered in fc e A M and it is different 
depending on the sublattice to which the reference cell k belongs. For a = 1,2,3,4 
we have 



iie{-l,0,+l}Cfc+ie2+ie 3 > k 6 A M , 

u j'e{-i,o,+i}Cfc+je 3 , k e A M , 

Cfc 7 k E A M , 

Uie{-i,o,+i}Cfc+ie 2 , k s A M . 



Figure [52] depicts the index sets A±{k) for the reference cell k belonging to A M for 
<> 1 1- 

General formulation. At this point we proceed with the general formulation for given 
a of the relevant quantities which are the reference measure v" a (da a ), the error term 

with k 6 Am and the sets A a (k) and B a (k), with the latter being the relevant 
boundary of A a . The index a indicates the sublattice we are considering. 

Definition 5.1. The sets A a (k) and B a (k) for k E A M are 

A a (k) = U/ : ||/_ fc || = i i p(;)< Q ,C; , B a (k) = Ui.||{_fc||=i, p (j)> a Cj . 



Definition 5.2. Given a = 1, . . . , 2 d we define the normalized Bernoulli measure 
on A M 

*% a (do»)= n (5-8) 

keAi, 

where 

= 7(A ( n.„> a . n(k v Z(A a (k)/{k}-a>^(k)) J] PM° a ). 



Z(A a (k):cr><*: V (k)) 



(5.9) 
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As we have seen in Step 3 we have two kinds of error terms in particular, 
those with k £ A M and others with k ^ A M . In order to describe the latter we need 
to introduce additional notation. 

For a = 1, . . . , 2 d we denote by T a the family of parallel hyperplanes of dimension 
d — 1 orthogonal to e a+ i passing through all the points k £ A M . Note that for any 
a, we have that Am = T Q U (r a + e Q +i). In the next definition we introduce a 
new parameter e a (k) £ {±1} depending on whether we should perform gluing or 
unfolding as discussed before. This is determined as follows: for fixed a £ Am let 
d(a, /3) be the distance between the sublattices A M and A M in the metric ||ck — /3||cx> = 
Si=i \ a i ~ Moreover, we can find orthogonal vectors {uj}j=i,...,d( a , / 8) and a family 
of signs {ej}j=i,...,d( a , / 8) such that 

d(a,0) 

A M = A M + 7(a. 0) with l(a, P) = ^ ej v,- . 
Note also that |7(a,/3)| = rf(a, f3). Then the exponents e a (fc) with p{k) = (3 are given 

by 

e a (k) := 

Furthermore, we denote by 7(0:, /?)) the affine hyperplane of codimension \ j(a 7 /3) | 
orthogonal to the connecting vectors {vj}j=i t ... l \-y( a ,p)\ an d passing through the point 
k 

Y(k,j(a,P))^n^ M Y(k, Vj ), 

where Y(k,v) is the hyperplane of dimension d — 1 passing through k and being 
perpendicular to the vector v. From the set of coarse-lattice points belonging to 
Y{k, v) we define the corresponding set by 



y(k, 7(a, /3)) := U ;e y (fc!7(Q!(3) )C; . 
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Then, letting k S A^ f , for I such that C; C dCk and with I e A^ /; for some /3, we 
define 

A a (O=( ifp(0>2^ 
[A Q (fc) n^G,7(a,/3)) otherwise, 

B a (0 = Ba(fe)ny(J,7(«./8))- (5-11) 

With the above definitions we can determine the error terms in the general expansion. 
Definition 5.3. For any k £ and for k €T a the error terms are given by 

$« = _! + Z (Mk); a >a -Mk))Z{A a+ i{k);a> a+1 - 1 n{k)) 
Moreover, if k G T a + e Q+ i and fc ^ A^/ 1 we /lave: 

1 -e a (k) 

>a+l. , '* 



Z(^ a (fc);(7 >Q ;7 ? (fc))Z(A a+ i(fc);a> a + 1 ; T? (fc)) 

( S 'fc- eQ + 1 ,e a+ i ^) ( S *+e„ + i ,e a+ i 



Furthermore, ifk € A^ 1 we replace Z(A*+i(fc); cr >Q!+1 ) 6j/ Z(A a+1 (k)/{k}; <7 >a+1 ). 

From Proposition 2.5.1 in [3 1 j we have that the general d-dimensional formulation 
of the a priori measure induced by the short-range interactions is 

e -*W Y[ P k (da)=R s (r ] )A(a)u(da\r ) ), 
fee A M 

where we have the following factors 

(i) a product of partition functions (depending only on the coarse-grained variable 

77) over finite sets of coarse cells with supports A 2 d(k), with k e A^ and 
<> 1 

RS M~Ii IT [z(A 2 4k); V(k)r^ k) ] , (5.12) 

(ii) error terms in the form of a gas of polymers (with the only interaction to be a 

hard-core exclusion) 

A(a):=n [] IT (! + **)> 

a=lj<2<*C«) fcgA^ 

(iii) a reference measure induced by only the short-range interactions once we neglect 

the reference system and the error terms 



v{da\rf) := v 2 



'>2">1 



With this expansion for the short range interactions, going back to the general 
strategy presented in Section [3J if we also consider the long-range contribution from 
\, we obtain 



e -/3H M (v) - / e -l3H l Ne -l3H 



k 

= e-^ m ^R( V ) [ e-^-^Aio-Mdaln) , 
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which implies that 

H M (v) =i/^°)(ry)-log J R(7 7 )-ilogE l/ [e-^«- S, ' (0) )A(a)|ry]. (5.13) 

5.1. Cluster expansion and effective interactions. The goal of this section 
is to expand the term ~R v [e~^^ HN ^~^ ' ° A{a)\rj\ in (|5.13[) into a convergent series 
using a cluster expansion. By the construction given previously the terms in A(a) 
are already in the form of a polymer gas with hard-core interactions only. For the 
long-range part we first write the difference H l N (a) — H l '^{rf) as 



H l N (a) - H l ^(r,) = ]T A kl J (a) , where 

k<l 

& u J{a);=-\ ]T (J(x-y)-J(k,l))a(x)a(y)(2-5 kl ). (5.14) 



2 .... 



We also define f k i(&) '■— e P^iJi^) _ i an( j we obtain 



_ (5-15) 

We define the polymer model which contains combined interactions originating 
from both the short and long-range potential. By expanding the products in (|5.15|) 
we obtain terms of the type 

p <? 

II II -^' m * where k ^ li ' m * 6 ^ M and G {1, . . . , 2 d } 

j=i i=i 

for some p and g. The factors are functions of the variables which are on the 
boundary of the corresponding sets A aj (fcj). This boundary is described by the set 

[B a+ i{k) it k e r a + e a+1 . 

Furthermore, since the measure v{da\rf) is not a product measure but instead a com- 
position of measures each one parametrized by variables which are integrated by the 
next measure, we need to create a "safety" corridor around the sets Cq depending 
on the level of a. This is given in the next definition. For a given integer /3, with 
1 < P < 2 d - a we define 

Cp(k) = U £l ,..., £/3e{±1} /3 U i:Ci cd(c fc+£1Co+1+ ... +e|3CQ+|3 ) >P (;)>a+/3 Ci , (5.17) 

Then for given a <E {1, . . . , 2 d } we call a "bond" of type C a the set 

C a (k) = Ufy*C${k). (5.18) 

With this definition, any factor has a region of dependence which is given by 
the bond C a (k). Similarly, for the factors ji i . mi originating from the long-range 
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interactions the initial domain of dependence is C/ i U C mi . However, due to the non- 
product structure of the measure we need to introduce a safety corridor in the same 
way. Given k £ Am for (3 an integer with 1 < (3 < 2 d — p{k) we define 

Cp(k) = U ei) .., ie?6{±1}( 3 Ul:C,Cd(C k+ t iea+1 + ... +tpea+l!l ),p(l)>P(k)+P Cl ■ ( 5 - 19 ) 

Then for a given fki we define 

C(k, I) = d^S^Cpik) u£i pW C (l) . (5.20) 
With a slight abuse of notation we define for Rq — {ki, . . . , k\R \} 

C(R ) = ufs* ki) Cp(ki) . (5.21) 

A bond Z will be either a Cj? bond for some a, k, called of type 1, or any subset 
of Am, we call it a bond of type 2. We say that two bonds l\ and I2 are connected 
if Zi n Z 2 7^ 0- We call a polymer R a set of bonds li,...,l p , Z p+1 where Z1, . . . , l p are 
bonds of type 1 and l p+ i is a bond of type 2, i.e., a possibly empty subset i?o C Am- 
A polymer is called connected if for all with l<i<j<p+l, there exists a 
chain of connected bonds in R joining Zj to Zj. For such a polymer i? we define its 
cardinality to be two integers, the first counting the number of bonds of the type 
1 and the second being the number of coarse cells included in the bond of type 2, 
i.e., card(i?) := (p, \Ro\). The support supp (R) of R is supp (R) = Zj where 
Z p+ i = C(Ro) (see (|5.2ip ). Let 1Z be the set of all such polymers. Two polymers 
Ri, i?2 are said to be compatible if Ri fl i?2 = and we write i?i ~ i?2- 

Given a polymer i? = Zi, . . . , l p , l p +\ we define the activity of R to be the function 
w : 1Z —> C given by 



iu(i?) = 




(5.22) 



where Qi is the collection of connected graphs on the vertices of Z (recall I C Am) and 
E(g) is the set of edges of the graph g. 

We define a new graph G on K which has the edge Ri-Rj if the polymers Ri 
and Rj are not compatible. We call G C 1Z completely disconnected if the subgraph 
induced by G on 67 has no edges. Let 

Vn = uJ2o{(fli> • • • >- R «) C ft : Vt 7^ j, Ri ~ iZ,- , } 
then the partition function Z can be written as 

z = E EH*)' 

which is the abstract form of a polymer model. Thus we can apply the general theorem 
of the cluster expansion once we check the convergence condition. The condition is 
stated as a theorem in [4]. 

Theorem 5.4 (j4|). Let a : 1Z -> R+. Consider the subset of C n 

Tk := {w(R), R e 1Z : VR e 1Z : \w(R)\e a( - R) < land 
2 (- iogC 1 - \ w (R')\e a{R ">)) < a(R)} ■ 
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Then on V^, \ogZ is well defined and analytic and 

\o g z= ci n mr) ir , 

IeX(K) flSsupp(I) 

where I = {Ir)r^-r, T(1Z) is the collection of all multi-indexes I, i.e., integer valued 
functions on 1Z, and 

_ 1 d lRl+ - +lR \*\logZ 

CI - I Rl \... I Rm ! Qlxi w{Rl) . . . Q** m w(R]nl) \W*)=°h ■ 

For the proof we refer to ,4 . Thus we need to check the condition of convergence. 
The following estimate for the long-range potential was proved in [35] . 

Lemma 5.5. Assume that J satisfies (B). Then there exists a constant C\ ~ 
q t llWHoo such that 

sup Y \&klJ(o-)\ <Ci, (5.23) 

for every a. 

For the short-range interaction we follow the analysis of [31] and we consider the 
following condition 

Condition 5.1. Let e be a vector in one of the directions of the lattice Am 
and Zu(A; er_, er + , r; r/v) be the partition function for the interaction U in the space 
domain A. We consider boundary conditions a± in the directions ±e andr in all other 
directions. Moreover, we impose multi- canonical constraints T}(k) for k G V C Aj\/ 
with A = UfegyCfc. For a given q > Tq, with \Ck\ = q d , the following inequality holds 



sup sup sup 

<y±,r A rjv 



Z u (A;a^ 7 a +7 T;i] V )Z u (A;0 7 0,T;riv) _ ^ 



Zu(A; 0, cr+, r; f] V )Zu{A; <r_, 0, r; n v ) 



< C 2 



where given the numbers r — 2 2rf [3(2 d+1 + l)] d , E = 2 d+1 + 1, and c > the upper 
bound C*2 satisfies 

rC 2 e cE < 1 . 



Notice that we work with the same condition as Condition Cl defined in [31 j . 
where in our notation L is q, yet similar analysis applies in order to prove convergence 
of the cluster expansion under the milder condition Condition C' L again as in [31 j - 
We skip the analysis of such issues since it goes beyond the goal of the present work. 
Furthermore, these conditions are related to the ones presented in [10] in order to 
ensure that a given system belongs to the class of completely analytical interactions. 
For further details we refer the reader to [3D] and [3J and to the references therein. 

Now we are ready to prove the convergence condition. 

Lemma 5.6. The set is nonempty. 
Proof: We take a(R) — c|supp (R)\, where c is a constant to be chosen later. Note 
that — log(l — x) < 2x, so it suffices to show that 

Y 2\w(R')\e a{R,) < a{R). 

R'ooR 
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Suppose that the generic polymer R' is given by R! = 1%, . . . , l p , l p +\, for some p > 0, 
where lj = C a j(k'j) for j = 1, . . . ,p and l p +i — R' 0) with \R' \ — n for some n > 0. 
For \w(R')\ we have 

\w(R')\< [v(da\ V )]l\t>%\-\ J2 II 

i=l 9ES R / {fc,i}G-B( S ) 

By the graph-tree inequality we have that for all <r, 77 and with \R' Q \ = n 

9£5r' {M}e£(g) r»eT« {M}eT° 

where from Lemma [531 we have that C\ ~ ^-j— 1| VV^Ijoo. We also let sup CT \AkiJ(cr)\ < 
Aki with Aki = g 2[l jjj || VV"||ool(ft ,/); |j-k|< h ■ Moreover, from Condition l5.ll we have 

v(da\ V )f[\$%\<(C 2 y. 



Then for the activity w(R') we obtain 



\w(R')\ </32"e" Cl £ II A w(<?2F 

r°er„° {fc,Z}ST0 



Thus to satisfy the sufficient condition for the convergence of the cluster expansion 
we first bound the sum J2r<*r 

by 



sup |supp(i?')|^E £ 

card <h')=(p, 



keR ' p>0n>0 n':SUpp (R')D{fc}, 



The fixed coarse cell Cfe may belong to one of the l/s for j = 1, . . . ,p or to R' . In 
the first case we estimate the sum over R' by 

1 



^ r , 0-1)! ^ (// - I.)! ' 

l 2,--,h k 2 ,...,k n 

and in the second by 

y - y -J—. 



Next, for every tree r° we have that 



sup y rj A fe()fc . ^cr 1 . 



1= i,j^ T 

&2 >• ■ ->fen 



We use the Cayley formula X!t <e7; 1 = n " 2 ano ^ the ^ ac ^ that the cardinality of the 
sum ^2 l i- 'i3{fe> can be bounded by r p , where r = 2 2(i [3(2 rf+1 + l)] d is an upper bound 
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for the maximum number of C a bonds that can pass through a point, as showed in 
[3Tj . Taking into account all the above we obtain 

J2 2\w(R')\e a W < |supp (R')\J2Mn)J2Mp), 

R'ooR n>\ p>l 

where 

Mn) = , 1 ^, n n - 2 (0C 1 ) n - 1 2 n e n ^e cn and A 2 (p) = , 1 ., r p C%e cEp , 
[n — 1)1 (p — iy. 

and E is the upper bound for the cardinality of any bond C a (k), i.e., 

sup sup \C a (k)\ < E = 2 d+1 + 1 . 

a=l,...,2 d feeA«j 

For C\ and C2 sufhciently small the two series converge to a finite number and we 
choose c to be this number. □ 
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